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Abstract 


The  SOC  computer  code,  designed  to 
calculate  the  effects  of  detonations  of  con¬ 
tained  explosions  on  the  surrounding 
medium, was  used  to  calculate  a  source 
function  called  the  reduced  displacement 
potential.  The  input  to  the  code  was  an 
experimentally  determined  equation  of 
state  for  each  of  the  following  frozen 
earth  materials  (at  -10°  C):  dry  and  ice- 
saturated  sand,  partially  and  fully  ice- 
saturated  glacial  till,  and  ice.  The 
reduced  displaced  potential  was  convolved 
with  the  response  of  the  Benioff 
seismometer  and  appropriate  attenuation 
operators.  The  first  half -cycle  displace¬ 


ment  amplitudes  at  300  and  500  km  were 
compared  to  those  obtained  from  particle 
velocity  measurements  for  nuclear  events 
in  four  mediums.  These  results,  scaled 
to  5  kt  and  normalized  to  tuff,  are: 
tuff  1,  granite  2.43,  salt  2.63,  and 
alluvium  0.18,  for  actual  nuclear  events; 
and  ice  0.78,  100%  saturated  sand  0.56, 
100%  saturated  till  0.67,  and  57%  saturated 
till  0.33,  for  simulated  nuclear  events  in 
these  mediums.  The  latter  value  for 
partially  saturated  till  assumed  a  com- 
pressional  velocity  of  2400  m/ sec  for  the 
medium.  A  lower  assumed  velocity  re¬ 
duces  the  values. 


Introduction 


Over  the  past  few  years,  considerable 
attention  has  been  focused  on  the  study 
of  seismic  signals  emanating  from  under¬ 
ground  nuclear  explosions.  The  primary 
interest  in  this  study  has  been  to  obtain 
an. understanding  of  the  attenuation  of  the 
signal  between  its  point  of  origin  and  re¬ 
cording  stations  tens  to  thousands  of 
kilometers  away. 

It  is  probably  safe  to  say  that  the 
seismic  signals  have  been  recorded  from 
every  underground  nuclear  detonation  that 
has  taken  place.  Analyses  of  these 
seismic  data  have  shown  that  the  amplitude 


of  the  signal  of  the  arrival  of  the  first 
cycle  and  a  half  is  dependent  on  the  energy 
released  and  the  propagation  distance  be¬ 
tween  source  and  station,  and  is  also 

strongly  dependent  on  the  shot  environ- 

1  2 
ment.  Werth,  Herbst,  and  Springer 

have  studied  the  factors  affecting  the 
first  arrival  of  seismic  waves  from  under¬ 
ground  nuclear  explosions  recorded  at 
distances  of  96  to  600  km.  At  these 
distances,  the  first  seismic  signal  to 
arrive  is  a  head  wave,  a  conical  wave 
arising  from  critical  refraction  at  the 
Mohorovicic  discontinuity.  Consequently, 
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Werth,  Herbst,  and  Springer  based  their 
technique  on  Zvolinskii's  head  wave 
theory,  which  relates  the  displacement 
amplitude  of  the  conical  wave  to  a  source 
function  called  the  reduced  displacement 
potential.  The  reduced  displacement 
potential  can  be  obtained  by  assuming 
that  the  motion  near  an  explosion,  but 
beyond  the  elastic  range  r  ,  can  be  rep¬ 
resented  by  the  elasticity  equation  for  a 
spherically  symmetric,  infinite  material. 
Then  a  "reduced"  displacement  potential 
A(t)  can  be  defined  by 

,  .v  _  9  A(t) 
u(re’ 1  "  8r  “7^  * 

where  r  =  t  -  r/c,  u(rg,  t)  is  the  displace¬ 
ment,  and  c  is  the  compressional  velocity 
of  the  medium.  The  distant  signal  is 
given  by  s(r)  =  6A(t),  where  9  is  some 
operator  representing  the  earth,  which 
corresponds  to  the  impulse  response, 
attenuation,  and  recording  instrument 
effects.5 

2 

Werth,  Herbst,  and  Springer  have 
shown  that  good  agreement  exists  between 
theory  ant.  measurement  for  several 
nuclear  events  in  tuff,  and  in  another  paper 
Werth  and  Herbst4  have  shown  that  when 
particle  displacement  data  just  beyond  the 
region  of  inelastic  medium  response  for 
nuclear  events  in  four  mediums  are  used 
to  calculate  the  reduced  displacement 
potentials,  the  resultant  amplitudes 
scaled  to  5  kt  and  normalized  to  tuff  are: 
tuff  1,  alluvium  0.25,  granite  1.11,  and 
salt  1.61 — demonstrating  the  effect  of  the 
shot  medium. 

Computer  codes  developed  at  Lawrence 
Radiation  Laboratory  have  been  used  to 
calculate,  with  considerable  success,  the 
effects  of  detonations  on  the  surrounding 


medium.5,6  Recently  the  codes  have  been 
improved  by  the  incorporation  of  a  more 
realistic  failure  model  and  better  mate¬ 
rial  descriptions. ^  With  adequate 
material  descriptions,  obtained  in  situ 
and  on  samples  recovered  from  the  shot 
environment,  excellent  agreement  between 
experimental  measurements  and  calcula¬ 
tions  is  consistently  obtained. 

The  intense  shock  generated  by  an 
underground  nuclear  explosion  vaporizes, 
fuses,  fractures,  and  permanently  dis¬ 
places  the  close-in  rock.  The  outgoing 
seismic  wave  is  essentially  spherical  and 
degenerates  from  a  strong  shock  wave  to 
an  elastic  wave.  The  equation  (in  finite 
difference  form)  relating  particle  dis¬ 
placement  time  history  at  some  radius  rg 
in  the  elastic  region  is  then  used  to  cal¬ 
culate  the  reduced  displacement  potential. 
The  reduced  displacement  potential  can 
be  calculated  for  any  complicated  u(t) 
from  outside  the  inelastic  region.  As 
long  as  no  characteristic  time  or  distances 
are  involved,  the  potential  formed  from 
the  displacement  history  of  an  explosion 
yield  (W^)  may  be  cube-root  scaled  to 

another  yield  (W9)  in  the  same  medium. 

,  ,.  ^  11 
The  scaling  goes  as 

A2(t)  =  A1(t)  (Wg/Wj) 

t2  .  Tj  (W^w//3 

The  continuing  program  in  the  seismic - 
decoupling  studies,  funded  in  part  by  the 
Advanced  Research  Project  Agency 
(ARPA),  includes  calculations  and  mate¬ 
rial  properties  studies.  Previous  studies 
suggest  that  materials  with  high  dry 
porosities  are  good  decoupling  mediums  in 
that  very  little  energy  is  transmitted  out 
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of  the  medium  into  the  seismic  wave. 
Other  materials  may  also  be  good  de¬ 
coupling  mediums  because  of  unique 
properties,  i.  e.,  phase  transitions  that 
take  place  under  shock  conditions. 


Perennially  frozen  ground,  or  perma¬ 
frost,  is  a  common  phenomenon  which 
underlies  a  major  portion  of  the  land 
areas  above  the  arctic  circle  to  depths  of 
as  much  as  500  m. 


$ 

,  Experimental 
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CHOICE  OF  MATERIALS 

The  existence  of  permafrost  is  deter¬ 
mined  by  thermal  conditions  and  is 
essentially  independent  of  the  composition 
of  earth  materials.  Virtually  any  type  of 
geological  material  may  be  perennially 
frozen.  The  properties  of  permafrost 
are  markedly  affected  by  the  presence  of 
ice  in  varying  quantaties,  which  indurates 
otherwise  unconsolidated  sediment.  There¬ 
fore,  the  natural  variability  of  earth  ma¬ 
terials  associated  with  varying  ice  volume 
makes  it  difficult  to  find  or  make-up  a 
typical  material.  The  samples  on  which 
tests  were  conducted  were  chosen  and 
made  up  by  the  U.  S.  Army  Cold  Regions 
Research  and  Engineering  Laboratory 
(CRUEL).12  The  samples  were  prepared 
at  varying  ice  contents  and  represent  a 
wide  range  of  naturally  occurring  conditions. 
The  sample  materials  are  (1)  Ottawa 
banding  sand,  a  well  rounded,  well  sorted, 
quartz  sand  with  a  median  diameter  of 
100  ju;  (2)  West  Lebanon  glacial  till,  a 
multiminer alic,  very  poorly  sorted, 
extremely  well  graded  material  with  no 
particle  larger  than  145  /u  and  a  median 
diameter  of  36  ju;  and  (3)  pure  polycrystal¬ 
line  and  single  crystal  ice,  the  matrix  ma¬ 
terial  in  permafrost. 


ISOTHERMAL  COMPRESSIBILITY 

The  technique  used  was  the  same  as 
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that  described  by  Stephens.  All  measure¬ 
ments  were  made  at  -10°C  (14°F). 

Figures  la-e  show  the  isothermal  com¬ 
pressibility  of  Ottawa  sand  with  varying 
degrees  of  saturation.  Figures  2a  and 
2b  show  the  data  for  glacial  till  at  57% 
and  100%  saturation.  Figure  3  shows  the 
compressibility  of  ice.  In  each  case 
both  the  loading  and  unloading  curves 
are  given. 

DYNAMIC  HUGONIOT  MEASUREMENTS 

Stanford  Research  Institute  (SRI)  made 
dynamic  Hugoniot  equation-of-state 
measurements  on  samples  with  the  same 
specifications  as  those  used  for  the  hydro¬ 
static  compressibility  measurements  of 
1 3 

CRREL.  Table  I  summarizes  the  re¬ 
sults.  All  tests  were  conducted  with  an 
initial  shot  temperature  of  -10°C. 

TRIAXLAL  STRENGTH  MEASUREMENTS 

Triaxial  and  unconfined  compressive 
strength  tests  were  conducted  on 
four  frozen  materials:  saturated  fine¬ 
grained  sand,  saturated  and  partially 
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Pressure  —  kbar  Pressure  —  kbar  Pressure  —  kbar 
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.  1.  Isothermal  compressibility  of  frozen  Ottawa  banding  sand  at  -10°C  for  five 
degrees  of  saturation  (S). 


Fig.  2.  Isothermal  compressibility  of  frozen  West  Lebanon  glacial  till  at  -10°  C  for 
58%  and  100%  saturation. 
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Table  I.  Dynamic  Hugoniot  equation-of-state  measurements  by  SRI  (-10°C). 


Material 

Dry  sand 

Pn  =  1.65  g/cm3 


Sand 

(20%  sat.) 

=  1.72  g/cm3 


Shock  velocity 
( cm/fusec ) 


Sand 

(50%  sat.)  c 
p  =  1.84  g/cm" 


0.407 

0.405 

0.531 

0.639 

0.298 

0.434 

0.424 

0.539 

0.599 

0.643 

0.340 

0.466 

0.532 

0.573 

0.637 


Particle  velocity 
(cm/jusec) 

0.217 

0.255 

0.313 

0.381 

0.114 

0.205 

0.212 

0.304 

0.349 

0.374 

0.111 

0.198 

0.279 

0.308 

0.344 


Pressure 

(kbar) 

Final 

volume 

(cm3/g) 

146 

0.2823 

150 

0.2097 

274 

0.2486 

404 

0.2449 

58 

0.3590 

155 

0.3071 

156 

0.2906 

282 

0.2527 

362 

0.2428 

413 

0.2431 

69.5 

0.3661 

170 

0.3127 

273 

0.2584 

325 

0.2509 

403 

0.2499 
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Table  I  (Continued) 


Final 

Shock  velocity 

Particle  velocity 

Pressure 

volume 

Material 

(cm/psec) 

(cm/psec) 

(kbar) 

(cm^/g) 

Sand 

0.377 

0.103 

76 

0.3708 

(100%  sat.) 
p  =  1.96  g/cm.3 

0.400 

0.101 

80 

0.3813 

0.386 

0.105 

80 

0.3714 

0.510 

0.201 

201 

0.3091 

0.577 

0.267 

303 

0.2736 

0.563 

0.271 

300 

0.2646 

0.596 

0.285 

339 

0.2632 

0.674 

0.324 

429 

0.2649 

0.677 

0.331 

440 

0.2603 

0.704 

0.352 

487 

0.2550 

West  Lebanon 

0.419 

0.208 

164 

0.2695 

glacial  till  (dry) 

p  =  1.86  g/cm.3 

0.417 

0.216 

167 

0.2589 

0.562 

0.294 

308 

0.2552 

0.657 

0.365 

451 

0.2393 

West  Lebanon 

0.359 

0.103 

76 

0.3470 

glacial  till 

0.352 

0.107 

77 

0.3347 

(57%  sat.) 

0.2958 

pQ  =  2.05  g/ cm^ 

0.482 

0.190 

187 

0.472 

0.196 

191 

0.2853 

0.570 

0.265 

309 

0.2610 

0.591 

0.263 

319 

0.2703 

0.605 

0.282 

354 

0.2593 

0.677 

0.319 

445 

0.2581 

0.658 

0.329 

444 

0.2440 

0.707 

0.345 

499 

0.2500 

West  Lebanon 

0.390 

0.100 

83 

0.3466 

glacial  till 
(100%  sat.) 

0.525 

0.194 

218 

0.2939 

0.2636 

pQ  =  2.145  g/cm* 

0.596 

0.259 

332 

0.701 

0.311 

467 

0.2593 

0.754 

0.334 

542 

0.2596 

Single -crystal  ice 

0.295 

0.130 

35.5 

0.6081 

po  =  0.917  g/crm* 

0.321 

0.125 

37.5 

0.666 

0.612 

0.319 

180 

0.5212 

0.755 

0.394 

274 

0.5207 

Polycrystalline  ice 

0.323 

0.125 

37.5 

0.6687 

pQ  =  0.917  g/cm^ 

0.500 

0.230 

106 

0.5889 

0.672 

0.336 

207 

0.5445 

0.781 

0.419 

301 

0.5055 

J 


v 
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saturated  glacial  till,  and  polycrystalline 
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ice.  All  tests  were  conducted  at  -10°C. 
For  frozen  soil  the  prevailing  mode  of 
failure  is  plastic  in  nature.  In  addition  to 
being  temperature  sensitive,  the  mode  of 
failure  is  strain  rate  sensitive — the 
frozen  soils  become  more  brittle  with 
increasing  strain  rate.  The  maximum 
value  of  mean  stress  obtainable  was 
limited  by  the  triaxial  apparatus  available. 

Since  the  matrix  material  in  frozen 
soils  is  ice,  strengths  are  strongly  in¬ 
fluenced  by  the  properties  of  the  ice 
matrix.  Figure  4  shows  the  results  of 
measurements  on  several  frozen  materi¬ 
als. 

SOC  CODE  CALCULATIONS 

In  the  measurement  program,  effec¬ 
tively  nine  different  materials  were 
tested, when  one  considers  that  each 
saturation  is  a  different  material.  For 
this  study,  five  materials  were  chosen 
as  representative  of  the  behavior  of 
permafrost.  These  are  (1)  dry  sand, 

(/> 


D 
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Fig.  4.  Triaxial  shear  strength  for  sev¬ 
eral  frozen  materials  at  -10°C. 


which  has  very  large  shear  strength  under 
high  average  stresses  and  no  strength 
at  zero  average  stress;  (2)  100%  saturated 
sand;  (3)  57%  saturated  glacial  till; 

(4)  100%  saturated  glacial  till;  and  (5)  the 
matrix  material,  ice.  The  calculations 
assumed  that  a  5-kt  nuclear  detonation 
had  taken  place  in  a  homogeneous  infinite 
medium  at  a  depth  of  burst  of  400  W1/3  ft 
(684  ft).  Inaddition,  a  5-kt  calculation  was 
made  for  ice  at  800  W"^3  depth  of  burst. 
The  equations  of  state  for  the  five  mate¬ 
rials  were  constructed  from  the  experi¬ 
mental  data  for  those  materials. 

Figures  5,  6,  7,  8  and  9  are  the  P-V 
curves  used  for  dry  sand,  100%  saturated 
sand,  57%  saturated  glacial  till,  100% 
saturated  glacial  till,  and  ice  respectively. 
The  lower  parts  of  the  curves,  up  to 
28  kbar,  are  identical  to  those  shown  in 
Figs.  1,  2,  and  3,  including  all  the  phase 
transitions  and  the  unloading  curves. 

Figures  lOa-e  show  the  shear  strength 
versus  mean  stress  relationships  used  in 
the  calculations.  The  points  are  identical 
to  those  shown  in  Fig.  4. 

Table  II  shows  various  values  used  in 
the  calculation.  pQ  is  the  bulk  density, 
k  is  the  elastic  bulk  modulus,  and  v0  is 
the  elastic  compressional  velocity  (m/sec). 
P1  and  P2  are  the  limits  of  mean  stress, 
where  it  was  assumed  the  materials 
behaved  bnttley;  above  P2  the  material 
was  assumed  to  be  plastic.  IC  ___  is  the 
maximum  allowable  shear  strength  after 
brittle  failure;  B  is  the  slope  of  the  maxi¬ 
mum  shear  strength  versus  mean  stress 
curve  after  brittle  failure,  where 
K  £  Bp  s  Kmax;  and  aelag  and  crfailed 
are  the  Poisson's  ratios  before  and  after 
failure  respectively. 
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Fig.  5.  Hugoniot  equation  of  state  for  dry  Fig.  6.  Hugoniot  equation  of  state  for  100% 

Ottawa  sand  at  -10°C.  saturated  Ottawa  sand  at  -10°C. 


Fig.  7.  Hugoniot  equation  of  state  for  Fig.  8.  Hugoniot  equation  of  state  for  100% 

57%  saturated  glacial  till  at  -10°C.  saturated  glacial  till  at  -10°C. 


Pressure  —  Mbar 


0.5 


3 

Volume — cm  / g 

Fig.  9.  Hugoniot  equation  of  state  for 
poly  crystalline  ice  at  -10°C. 

The  values  of  bulk  modulus  k  and 
elastic  velocity  vg  were  determined  either 
from  handbook  values  (dry  sand  and  ice) 
or  from  the  initial  slope  of  the  isothermal 
compressibility  curve,  using  the  elastic 
equations 

<ve)2po  =  k  +  f^ 


where  /u  is  the  elastic  shear  modulus,  and 


For  the  57%  saturated  till,  a  velocity  of 
350  m/sec  is  obtained  from  the  iso¬ 
thermal  compressibility  curve  to  125  bars, 
with  an  assumed  Poisson's  ratio  of  0.35. 
The  fact  that  the  100%  saturated  till  gives 
a  velocity  of  3500  m/ sec,  makes  it  dif¬ 
ficult  to  believe  that  a  factor-of-10 
difference  in  velocity  occurs  between 
these  two  saturations;  consequently  an 
additional  calculation  was  made  using  a 
reported  value  of  2400  ra./sec}^ 

RESULTS  OF  CALCULATIONS 

Table  III  lists  some  of  the  results  of 
the  six  calculations  made  on  five  mate¬ 
rials.  Two  calculations  were  made  on 
ice  to  test  the  effect  of  depth  of  burst  on 
the  various  calculated  quantities. 

The  results  indicate  that,  in  general, 
the  static  state  reduced  displacement 
potential,  is  related  directly  to  the  cavity 
radius.  The  steady  state  reduced  dis¬ 
placement  potential  is  not  given  for  dry 
sand.  Dry  sand  was  not  allowed  to  under¬ 
go  brittle  failure,  and  continues  to  yield 


Table  II.  Parameters  used  in  SOC  code  calculations. 


po 

(g/  cm3) 

k 

(kbar)  ( 

ve  ^1  -^2  ^max 

m/ sec)  (kbar)  (kbar)  (kbar) 

B 

aelas 

°failed 

A 

Dry  sand 

1.65 

7.09 

833 

0 

0 

6 

0.5 

0.3 

0.3 

100%  saturated 
sand 

1.985 

107.09 

2800 

0 

0.05 

0.2 

0.75 

0.35 

0.35 

1 

57%  saturated  till 

2.03 

1.72 

80.95 

350 

2400 

0 

0.05 

0.16 

0.7 

0.35 

0.35 

Ice 

0.917 

3600 

0 

0.2 

0.16 

0.4 

0.35 

0.35 
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Fig.  10a.  Triaxial  shear  strength  versus  mean  stress  for  dry  Ottawa  sand  at  -10  C. 


Fig.  10b.  Triaxial  shear  strength  versus  mean  stress  for  100%  saturated  sand  at  -10  C 


plastically  to  zero  mean  stress.  Conse¬ 
quently  the  calculated  reduced  displace¬ 
ment  potential  for  dry  sand  continues  to 
decrease  with  distance.  The  reduced 
displacement  potential  is  a  strong  function 
of  the  porosity  of  the  medium.  For  a 
partially  saturated  material  with  dry 


porosity  like  the  till,  the  reduced  displace¬ 
ment  potential  is  not  meaningful  except 
for  that  particular  saturation. 

Figure  11  shows  the  calculated  peak 
particle  velocity  versus  distance  for  the 
five  materials  at  overburden  (pgh)  pres¬ 
sures  of  the  208-m  depth  of  burst. 


jr 
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Fig.  lOe.  Triaxial  shear  strength  versus  mean  stress  for  poly  crystalline  ice  at  -10°C. 


Fig.  11.  Calculated  peak  particle  velocity 
as  a  function  of  distance  for  the 
five  materials  studied  (DOB  = 

208  m). 

Figure  12  shows  the  peak  radial  stress  as 
a  function  of  distance  for  the  same  calcu¬ 
lations.  Figure  13  shows  the  radial  stress 
pulse  shape  at  30  m  for  each  calculation. 
Figure  14  shows  the  reduced  displacement 
potential  for  those  calculations  reported 
in  Table  III  at  208 -m  depth  of  burst.  The 


Fig.  12.  Calculated  peak  radial  stress  as 
a  function  of  distance  for  the  five 
materials  studied  (with  over¬ 
burden  subtracted  out)  (DOB  = 
208  m). 

calculation  made  with  ice  at  416 -m  depth 
of  burst  shows  the  effect  of  depth  of  burst 
on  the  reduced  displacement  potential. 

ANALYSIS  OF  RESULTS 

The  analysis  made  is  based  on  a 
comparison  of  the  calculated  first  half- 
cycle  displacement  as  reported  by  Werth 

4 

and  Herbst.  They  calculated  reduced 
displacement  potential  from  particle 


/ 


"X 

) 


Ti  me  —  msec 


Fig.  13.  Calculated  radial  stress  pulse  shapes  for  5-kt  detonations  in  four  frozen 
earth  materials  at  30  m  from  working  point  (DOB  =  208  m). 

velocity  measurements  just  beyond  the 
region  of  plastic  deformations  for  nuclear 
explosions  in  four  mediums.  The  results 
were  scaled  to  5  kt  and  the  potentials  were 
convolved  with  the  response  of  the  Benioff 
instrument  with  an  appropriate  attenuation 
operator.  The  models  of  the  earth  were 
developed  for  each  explosion,  and  the  re¬ 
sults  were  scaled  to  tuff  as  shown  in 
Table  IV. 

16 

Roberts  calculated  first  half- cycle 
displacements  versus  distance  for  each 
reduced  displacement  potential  shown  in 
Fig.  13.  He  used  the  values  of  density, 
and  compressional  and  shear  velocity 
shown  in  Table  III,  and  the  crustal  model 
for  the  Nevada  Test  Site.  Table  V  gives 
the  results  of  these  calculations. 

With  these  values,  a  direct  compari¬ 
son  with  the  four  nuclear  explosive 


Time  — sec 


Fig.  14.  Calculated  reduced  displacement 
potentials  for  four  different 
media. 
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Table  III.  Results  of  SOC  code  calculations. 


Material 

Final 

cavity 

radius 

(m) 

Radius  of 
brittle 
failure 
(m) 

Steady 

state 

RDPa 

(m3) 

Overburden 

pressure 

(bars) 

Equivalent 
depth  of 
burst 
(m) 

Dry  sand 

23.96 

— 

— 

34 

208  (684  ft) 

100%  sat.  sand 

20.12 

156.7 

1880 

40 

208 

57%  sat.  till 

29.88 

76.1 

— 

41 

208 

100%  sat.  till 

19.41 

252.5 

1680 

44 

208 

Ice 

26.53 

268.2 

4300 

20 

208 

Ice 

23.95 

206.0 

3100 

40 

416  (1368  ft) 

ct 

Reduced  displacement  potential. 


Table  IV.  Five-kiloton  source  with  attenuation  and  crustal  models  (after  Werth  and 
Herbst).4 


Tuff 

Alluvium 

Granite 

Salt 

First  half -cycle  amplitude 
at  300  km 

92  m fj, 

16  m fj, 

224  mp 

242  m/u 

Ratio  to  tuff 

1 

0.17 

2.43 

2.63 

First  half-cycle  amplitude 
at  500  km 

19  m/i 

3.5  nyj 

43  m^u 

40  m ifj, 

Ratio  to  tuff 

1 

0.18 

2.26 

2.11 

Table  V.  Five-kiloton  source  with  attenuation  and  crustal  models. 


_ First  half -cycle  amplitude  (mg) _ 

Glacial  till _ Sand  _ Ice 


Distance 

(km) 

57%  sat. 
ve  =  350 
m/  sec 

57%  sat. 
ve  =  2400 
m/sec 

100% 

sat. 

100% 

sat. 

pgh  =  40 
bars 

pgh  =  20 
bars 

200 

15.9 

95.7 

197.5 

163.8 

199.6 

2218 

300 

5.1 

30.7 

61.5 

51.4 

62.6 

71.6 

400 

2.2 

13.2 

25.7 

21.6 

26.9 

30.8 

500 

1.1 

6.8 

13.0 

11.0 

13.6 

15.9 

600 

0.6 

3.8 

7.2 

6.1 

7.7 

9.0 

700 

0.4 

2.4 

4.4 

3.7 

4.7 

5.6 

mediums  used  by  Werth  and  Herbst  was 
made.  Table  VI  shows  the  results. 

The  above  results  seem  to  indicate 
that  saturated  frozen  earth  materials  are 
poorer  decoupling  mediums  than  Nevada 


Test  Site  valley  alluvium.  An  interesting 
result  is  that  the  assumed  low  velocity 
for  partially  saturated  glacial  till 
produces  appreciably  lower  amplitudes. 
This  is  due  to  low  compressional  wave 


J 


} 


1 

) 


Table  VI.  Resultant  amplitudes  from  Table  V  normalized  to  tuff. 


Amplitudes  (mju) 

Glacial  till 

Sand 

Ice 

Tuff 

57%  sat. 
ve  =  350 
m/  sec 

57%  sat. 
ve  =  2400 
m/  sec 

100% 

sat. 

100% 

sat. 

pgh  =  40 
bars 

pgh  =  20 
bars 

First  half-cycle 
amplitude  at 
300  km 

92 

5.1 

30.7 

61.6 

51.4 

62.6 

71.6 

Ratio  to  tuff 

1 

0.06 

0.33 

0.67 

0.56 

0.68 

0.78 

First  half-cycle 
amplitude  at 
500  km 

19 

1.1 

6.8 

13.0 

11.0 

13.6 

15.9 

Ratio  to  tuff 

1 

0.06 

0.36 

0.68 

0.58 

0.72 

0.84 

velocity,  which  leads  to  an  unusually 
low  coupling  of  the  seismic  energy  out 
of  the  detonation  medium.  Another 
interesting  result  shown  in  Table  VI 
shows  the  effect  of  depth  of  burst  on 
the  predicted  amplitude. 

With  regard  to  the  possible  low 
coupling  of  partially  saturated  glacial 
till,  it  is  very  unlikely  that  partially 
saturated  material  will  be  found  at 
containment  depths  because  of  the 
plastic- viscous  nature  of  the  matrix 
material  (ice).  Because  of  its  ductile 
nature,  loose  dry  sand  could  prove  to 
be  an  excellent  decoupling  medium, 
provided  it  could  be  found  at  a  suffi¬ 
cient  depth  and  thickness. 

Figure  15  shows  a  comparison  of 
calculated  first  half-cycle  amplitudes 
(versus  distance)  of  the  frozen  materials 

4 

with  those  reported  by  Werth  and  Herbst. 

Fig.  15.  Comparison  of  calculated  first  half - 
^  cycle  displacement  versus  distance 

for  several  real  materials  and  sev¬ 
eral  hypothetical  materials.  Displace 
ments  are  normalized  to  5 -kt.  (Salt 
—  Gnome  Event;  granite — Hardhat 
Event;  tuff — Rainier  Event;  allu¬ 
vium —  Fisher  Event.) 
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